Transcription-replication interactions reveal bacterial genome regulation.

Organisms determine the transcription rates of thousands of genes through a few modes of regulation that recur across the genome1. In bacteria, the relationship between the regulatory architecture of a gene and its expression is well understood for individual model gene circuits2,3. However, a broader perspective of these dynamics at the genome scale is lacking, in part because bacterial transcriptomics has hitherto captured only a static snapshot of expression averaged across millions of cells4. As a result, the full diversity of gene expression dynamics and their relation to regulatory architecture remains unknown. Here we present a novel genome-wide classification of regulatory modes based on the transcriptional response of each gene to its own replication, which we term the transcription-replication interaction profile (TRIP). Analysing single-bacterium RNA-sequencing data, we found that the response to the universal perturbation of chromosomal replication integrates biological regulatory factors with biophysical molecular events on the chromosome to reveal the local regulatory context of a gene. Whereas the TRIPs of many genes conform to a gene dosage-dependent pattern, others diverge in distinct ways, and this is shaped by factors such as intra-operon position and repression state. By revealing the underlying mechanistic drivers of gene expression heterogeneity, this work provides a quantitative, biophysical framework for modelling replication-dependent expression dynamics.

Association of alpha globin gene copy number with exhaled nitric oxide in a cross-sectional study of healthy Black adults.

INTRODUCTION: The genetic determinants of fractional exhalation of nitric oxide (FeNO), a marker of lung inflammation, are understudied in Black individuals. Alpha globin (HBA) restricts nitric oxide signalling in arterial endothelial cells via interactions with nitric oxide synthase; however, its role in regulating the release of NO from respiratory epithelium is less well understood. We hypothesised that an HBA gene deletion, common among Black individuals, would be associated with higher FeNO. METHODS: Healthy Black adults were enrolled at four study sites in North Carolina from 2005 to 2008. FeNO was measured in triplicate using a nitric oxide analyzer. The -3.7 kb HBA gene deletion was genotyped using droplet digital PCR on genomic DNA. The association of FeNO with HBA copy number was evaluated using multivariable linear regression employing a linear effect of HBA copy number and adjusting for age, sex and serum immunoglobulin-E levels. Post-hoc analysis employing a recessive mode of inheritance was performed. RESULTS: 895 individuals were in enrolled in the study and 720 consented for future genetic research; 643 had complete data and were included in this analysis. Median (25th, 75th) FeNO was 20 (13, 31) ppb. HBA genotypes were: 30 (4.7%) -a/-a, 197 (30.6%) -a/aa, 405 (63%) aa/aa and 8 (1.2%) aa/aaa. Subjects were 35% male with median age 20 (19, 22) years. Multivariable linear regression analysis revealed no association between FeNO and HBA copy number (ß=-0.005 (95% CI -0.042 to 0.033), p=0.81). In the post-hoc sensitivity analysis, homozygosity for the HBA gene deletion was associated with higher FeNO (ß=0.107 (95% CI 0.003 to 0.212); p=0.045). CONCLUSION: We found no association between HBA copy number and FeNO using a prespecified additive genetic model. However, a post hoc recessive genetic model found FeNO to be higher among subjects homozygous for the HBA deletion.

Alpha globin gene copy number and hypertension risk among Black Americans.

BACKGROUND: Alpha globin is expressed in the endothelial cells of human resistance arteries where it binds to endothelial nitric oxide synthase and limits release of the vasodilator nitric oxide. Genomic deletion of the alpha globin gene (HBA) is common among Black Americans and could lead to increased endothelial nitric oxide signaling and reduced risk of hypertension. METHODS: Community-dwelling US adults aged 45 years or older were enrolled and examined from 2003 to 2007, followed by telephone every 6 months, and reexamined from 2013 to 2016. At both visits, trained personnel performed standardized, in-home blood pressure measurements and pill bottle review. Prevalent hypertension was defined as systolic blood pressure ≥ 140mmHg or diastolic blood pressure ≥ 90mmHg or anti-hypertensive medication use. Droplet digital PCR was used to determine HBA copy number. The associations of HBA copy number with prevalent hypertension, resistant hypertension, and incident hypertension were estimated using multivariable regression. RESULTS: Among 9,684 Black participants, 7,439 (77%) had hypertension at baseline and 1,044 of those had treatment-resistant hypertension. 1,000 participants were not hypertensive at baseline and participated in a follow up visit; 517 (52%) developed hypertension over median 9.2 years follow-up. Increased HBA copy number was not associated with prevalent hypertension (PR = 1.00; 95%CI 0.98,1.02), resistant hypertension (PR = 0.95; 95%CI 0.86,1.05), or incident hypertension (RR = 0.96; 95%CI 0.86,1.07). CONCLUSIONS: There were no associations between increased HBA copy number and risk of hypertension. These findings suggest that variation in alpha globin gene copy number does not modify the risk of hypertension among Black American adults.

Why Females Do Better: The X Chromosomal TLR7 Gene-Dose Effect in COVID-19.

A male sex bias has emerged in the COVID-19 pandemic, fitting to the sex-biased pattern in other viral infections. Males are 2.84 times more often admitted to the ICU and mortality is 1.39 times higher as a result of COVID-19. Various factors play a role in this, and novel studies suggest that the gene-dose of Toll-Like Receptor (TLR) 7 could contribute to the sex-skewed severity. TLR7 is one of the crucial pattern recognition receptors for SARS-CoV-2 ssRNA and the gene-dose effect is caused by X chromosome inactivation (XCI) escape. Female immune cells with TLR7 XCI escape have biallelic TLR7 expression and produce more type 1 interferon (IFN) upon TLR7 stimulation. In COVID-19, TLR7 in plasmacytoid dendritic cells is one of the pattern recognition receptors responsible for IFN production and a delayed IFN response has been associated with immunopathogenesis and mortality. Here, we provide a hypothesis that females may be protected to some extend against severe COVID-19, due to the biallelic TLR7 expression, allowing them to mount a stronger and more protective IFN response early after infection. Studies exploring COVID-19 treatment via the TLR7-mediated IFN pathway should consider this sex difference. Various factors such as age, sex hormones and escape modulation remain to be investigated concerning the TLR7 gene-dose effect.

An optimal growth law for RNA composition and its partial implementation through ribosomal and tRNA gene locations in bacterial genomes.

The distribution of cellular resources across bacterial proteins has been quantified through phenomenological growth laws. Here, we describe a complementary bacterial growth law for RNA composition, emerging from optimal cellular resource allocation into ribosomes and ternary complexes. The predicted decline of the tRNA/rRNA ratio with growth rate agrees quantitatively with experimental data. Its regulation appears to be implemented in part through chromosomal localization, as rRNA genes are typically closer to the origin of replication than tRNA genes and thus have increasingly higher gene dosage at faster growth. At the highest growth rates in E. coli, the tRNA/rRNA gene dosage ratio based on chromosomal positions is almost identical to the observed and theoretically optimal tRNA/rRNA expression ratio, indicating that the chromosomal arrangement has evolved to favor maximal transcription of both types of genes at this condition.

High-Throughput and Accurate Determination of Transgene Copy Number and Zygosity in Transgenic Maize: From DNA Extraction to Data Analysis.

It is vital to develop high-throughput methods to determine transgene copy numbers initially and zygosity during subsequent breeding. In this study, the target sequence of the previously reported endogenous reference gene hmg was analyzed using 633 maize inbred lines, and two SNPs were observed. These SNPs significantly increased the PCR efficiency, while the newly developed hmg gene assay (hmg-taq-F2/R2) excluding these SNPs reduced the efficiency into normal ranges. The TaqMan amplification efficiency of bar and hmg with newly developed primers was calculated as 0.993 and 1.000, respectively. The inter-assay coefficient of variation (CV) values for the bar and hmg genes varied from 1.18 to 2.94%. The copy numbers of the transgene bar using new TaqMan assays were identical to those using dPCR. Significantly, the precision of one repetition reached 96.7% of that of three repetitions of single-copy plants analyzed by simple random sampling, and the actual accuracy reached 95.8%, confirmed by T1 and T2 progeny. With the high-throughput DNA extraction and automated data analysis procedures developed in this study, nearly 2700 samples could be analyzed within eight hours by two persons. The combined results suggested that the new hmg gene assay developed here could be a universal maize reference gene system, and the new assay has high throughput and high accuracy for large-scale screening of maize varieties around the world.

Circulating tumor DNA-guided treatment with pertuzumab plus trastuzumab for HER2-amplified metastatic colorectal cancer: a phase 2 trial.

The applicability of circulating tumor DNA (ctDNA) genotyping to inform enrollment of patients with cancer in clinical trials has not been established. We conducted a phase 2 trial to evaluate the efficacy of pertuzumab plus trastuzumab for metastatic colorectal cancer (mCRC), with human epidermal growth factor receptor 2 (HER2) amplification prospectively confirmed by tumor tissue or ctDNA analysis ( UMIN000027887 ). HER2 amplification was confirmed in tissue and/or ctDNA in 30 patients with mCRC. The study met the primary endpoint with a confirmed objective response rate of 30% in 27 tissue-positive patients and 28% in 25 ctDNA-positive patients, as compared to an objective response rate of 0% in a matched real-world reference population treated with standard-of-care salvage therapy. Post hoc exploratory analyses revealed that baseline ctDNA genotyping of HER2 copy number and concurrent oncogenic alterations adjusted for tumor fraction stratified patients according to efficacy with similar accuracy to tissue genotyping. Decreased ctDNA fraction 3 weeks after treatment initiation associated with therapeutic response. Pertuzumab plus trastuzumab showed similar efficacy in patients with mCRC with HER2 amplification in tissue or ctDNA, showing that ctDNA genotyping can identify patients who benefit from dual-HER2 blockade as well as monitor treatment response. These findings warrant further use of ctDNA genotyping in clinical trials for HER2-amplified mCRC, which might especially benefit patients in first-line treatment.

Transcriptomic and rRNA:rDNA Signatures of Environmental versus Enteric Enterococcus faecalis Isolates under Oligotrophic Freshwater Conditions.

The use of enterococci as a fecal indicator bacterial group for public health risk assessment has been brought into question by recent studies showing that "naturalized" populations of Enterococcus faecalis exist in the extraenteric environment. The extent to which these naturalized E. faecalis organisms can confound water quality monitoring is unclear. To determine if strains isolated from different habitats display different survival strategies and responses, we compared the decay patterns of three E. faecalis isolates from the natural environment (environmental strains) against three human gut isolates (enteric strains) in laboratory mesocosms that simulate an oligotrophic, aerobic freshwater environment. Our results showed similar overall decay rates between enteric and environmental isolates based on viable plate and quantitative PCR (qPCR) counts. However, the enteric isolates exhibited a spike in copy number ratios of 16S rRNA gene transcripts to 16S rRNA gene DNA copies (rRNA:rDNA ratios) between days 1 and 3 of the mesocosm incubations that was not observed in environmental isolates, which could indicate a different stress response. Nevertheless, there was no strong evidence of differential gene expression between environmental and enteric isolates related to habitat adaptation in the accompanying mesocosm metatranscriptomes. Overall, our results provide novel information on how rRNA levels may vary over different growth conditions (e.g., standard lab versus oligotrophic) for this important indicator bacteria. We also observed some evidence for habitat adaptation in E. faecalis; however, this adaptation may not be substantial or consistent enough for integration in water quality monitoring. IMPORTANCE Enterococci are commonly used worldwide to monitor environmental fecal contamination and public health risk for waterborne diseases. However, closely related enterococci strains adapted to living in the extraenteric environment may represent a lower public health risk and confound water quality estimates. We developed an rRNA:rDNA viability assay for E. faecalis (a predominant species within this fecal group) and tested it against both enteric and environmental isolates in freshwater mesocosms to assess whether this approach can serve as a more sensitive water quality monitoring tool. We were unable to reliably distinguish the different isolate types using this assay under the conditions tested; thus, environmental strains should continue to be counted during routine water monitoring. However, this assay could be useful for distinguishing more recent (i.e., higher-risk) fecal pollution because rRNA levels significantly decreased after 1 week in all isolates.

Effect of EPSPS gene copy number and glyphosate selection on fitness of glyphosate-resistant Bassia scoparia in the field.

The widespread evolution of glyphosate-resistant (GR) Bassia scoparia in the U.S. Great Plains poses a serious threat to the long-term sustainability of GR sugar beet. Glyphosate resistance in B. scoparia is due to an increase in the EPSPS (5-enolpyruvyl-shikimate-3-phosphate) gene copy number. The variation in EPSPS gene copies among individuals from within a single GR B. scoparia population indicated a differential response to glyphosate selection. With the continued use of glyphosate in GR sugar beet, the effect of increasing glyphosate rates (applied as single or sequential applications) on the fitness of GR B. scoparia individuals with variable EPSPS gene copies was tested under field conditions. The variation in EPSPS gene copy number and total glyphosate rate (single or sequential applications) did not influence any of the reproductive traits of GR B. scoparia, except seed production. Sequential applications of glyphosate with a total rate of 2214 g ae ha-1 or higher prevented seed production in B. scoparia plants with 2-4 (low levels of resistance) and 5-6 (moderate levels of resistance) EPSPS gene copies. Timely sequential applications of glyphosate (full recommended rates) can potentially slow down the evolution of GR B. scoparia with low to moderate levels of resistance (2-6 EPSPS gene copies), but any survivors (highly-resistant individuals with ≥ 8 EPSPS gene copies) need to be mechanically removed before flowering from GR sugar beet fields. This research warrants the need to adopt ecologically based, multi-tactic strategies to reduce exposure of B. scoparia to glyphosate in GR sugar beet.

Gene copy-number changes and chromosomal instability induced by aneuploidy confer resistance to chemotherapy.

Mitotic errors lead to aneuploidy, a condition of karyotype imbalance, frequently found in cancer cells. Alterations in chromosome copy number induce a wide variety of cellular stresses, including genome instability. Here, we show that cancer cells might exploit aneuploidy-induced genome instability and the resulting gene copy-number changes to survive under conditions of selective pressure, such as chemotherapy. Resistance to chemotherapeutic drugs was dictated by the acquisition of recurrent karyotypes, indicating that gene dosage might play a role in driving chemoresistance. Thus, our study establishes a causal link between aneuploidy-driven changes in gene copy number and chemoresistance and might explain why some chemotherapies fail to succeed.

Genomic landscape and evolution of arm aneuploidy in lung adenocarcinoma.

For lung adenocarcinoma, arm aneuploidy landscape among primary and metastatic sites, and among different driver and frequently mutated gene groups have not been previously studied. We collected the largest cohort of LUAD patients (n=3533) to date and analyzed the profiles of chromosome arm aneuploidy (CAA), and its association with different metastatic sites and mutated gene groups. Our results showed distant metastasis (bone, brain, liver) were characterized by high CAA burden and biased towards arm losses compared to regional metastasis (pleura, chest) and primary tumors. Moreover, EGFR, MET, PIK3CA, PKHD1 and RB1 mutant groups were found to have high CAA burden, while those with BRAF, ERBB2 and KRAS mutations belonged to the low CAA burden group. Comparing EGFR L858R and EGFR 19del mutants, distinct CAA co-occurrences were observed. Network-based stratification with population based genomic evolution analysis revealed two distinct subtypes of LUAD with different CAA signatures and unique CAA order of acquisition. In summary, our study presented a comprehensive characterization of arm aneuploidy landscape and evolutionary trajectories in lung adenocarcinoma, which could provide basis for both biological and clinical investigations in the future.

Heuristics for Genome Rearrangement Distance With Replicated Genes.

In comparative genomics, one goal is to find similarities between genomes of different organisms. Comparisons using genome features like genes, gene order, and regulatory sequences are carried out with this purpose in mind. Genome rearrangements are mutational events that affect large extensions of the genome. They are responsible for creating extant species with conserved genes in different positions across genomes. Close species - from an evolutionary point of view - tend to have the same set of genes or share most of them. When we consider gene order to compare two genomes, it is possible to use a parsimony criterion to estimate how close the species are. We are interested in the shortest sequence of genome rearrangements capable of transforming one genome into the other, which is named rearrangement distance. Reversal is one of the most studied genome rearrangements events. This event acts in a segment of the genome, inverting the position and the orientation of genes in it. Transposition is another widely studied event. This event swaps the position of two consecutive segments of the genome. When the genome has no gene repetition, a common approach is to map it as a permutation such that each element represents a conserved block. When genomes have replicated genes, this mapping is usually performed using strings. The number of replicas depends on the organisms being compared, but in many scenarios, it tends to be small. In this work, we study the rearrangement distance between genomes with replicated genes considering that the orientation of genes is unknown. We present four heuristics for the problem of genome rearrangement distance with replicated genes. We carry out experiments considering the exclusive use of the reversals or transpositions events, as well as the version in which both events are allowed. We developed a database of simulated genomes and compared our results with other algorithms from the literature. The experiments showed that our heuristics with more sophisticated rules presented a better performance than the known algorithms to estimate the evolutionary distance between genomes with replicated genes. In order to validate the application of our algorithms in real data, we construct a phylogenetic tree based on the distance provided by our algorithm and compare it with a know tree from the literature.

Polyploid mitosis and depolyploidization promote chromosomal instability and tumor progression in a Notch-induced tumor model.

Ploidy variation is a cancer hallmark and is frequently associated with poor prognosis in high-grade cancers. Using a Drosophila solid-tumor model where oncogenic Notch drives tumorigenesis in a transition-zone microenvironment in the salivary gland imaginal ring, we find that the tumor-initiating cells normally undergo endoreplication to become polyploid. Upregulation of Notch signaling, however, induces these polyploid transition-zone cells to re-enter mitosis and undergo tumorigenesis. Growth and progression of the transition-zone tumor are fueled by a combination of polyploid mitosis, endoreplication, and depolyploidization. Both polyploid mitosis and depolyploidization are error prone, resulting in chromosomal copy-number variation and polyaneuploidy. Comparative RNA-seq and epistasis analysis reveal that the DNA-damage response genes, also active during meiosis, are upregulated in these tumors and are required for the ploidy-reduction division. Together, these findings suggest that polyploidy and associated cell-cycle variants are critical for increased tumor-cell heterogeneity and genome instability during cancer progression.

The oncological relevance of fragile sites in cancer.

Recent developments in sequencing the cancer genome have provided the first in-depth mapping of structural variants (SV) across 38 tumour types. Sixteen signatures of structural variants have been proposed which broadly characterise the variation seen across cancer types. One signature shows increased duplications and deletions at fragile sites, with little association with the typical DNA repair defects. We discuss how, for many of these fragile sites, the clinical impacts are yet to be explored. One example is NAALADL2, one of the most frequently altered fragile sites in the cancer genome. The copy-number variations (CNVs) which occur at fragile sites, such as NAALADL2, may span many genes without typical DNA repair defects and could have a large impact on cell signalling.

The Flora Compositions of Nitrogen-Fixing Bacteria and the Differential Expression of nifH Gene in Pennisetum giganteum z.x.lin Roots.

The flora compositions of nitrogen-fixing bacteria in roots of Pennisetum giganteum z.x.lin at different growth stages and the expression and copy number of nitrogen-fixing gene nifH were studied by Illumina Miseq second-generation sequencing technology and qRT-PCR. The results showed that there were more than 40,000~50,000 effective sequences in 5 samples from the roots of P. giganteum, with Proteobacteria and Cyanobacteria as the dominant nitrogen-fixing bacteria based on the OTU species annotations for each sample and Bradyrhizobium as the core bacterial genera. The relative expression and quantitative change of nifH gene in roots of P. giganteum at different growth stages were consistent with the changes in the flora compositions of nitrogen-fixing microbia. Both revealed a changing trend with an initial increase and a sequential decrease, as well as changing order as jointing stage>maturation stage>tillering stage>seedling stage>dying stage. The relative expression and copy number of nifH gene were different in different growth stages, and the difference among groups basically reached a significant level (p < 0.05). The relative expression and copy number of nifH gene at the jointing stage were the highest, and the 2-△△CT value was 4.43 folds higher than that at the seedling stage, with a copy number of 1.32 × 107/g. While at the dying stage, it was the lowest, and the 2-△△CT value was 0.67 folds, with a copy number of 0.31 × 107/g.

Evidence for frequent concurrent DCUN1D1, FGFR1, BCL9 gene copy number amplification in squamous cell lung cancer.

Non-small cell lung cancer (NSCLC) targeted therapies are mostly based on activating mutations and rearrangements which are rare events in Lung Squamous Cell Carcinomas (LUSC). Recently advances in immunotherapy have improved the therapeutic repository for LUSC, but there is still an urgent need for novel targets and biomarkers. We examined 73 cases of LUSC for relative copy number amplification of DCUN1D1, BCL9, FGFR1 and ERBB2 genes and searched for correlations with molecular alterations and clinicopathological characteristics. In our cohort BCL9 gene was amplified in 57.5 % of the cases, followed by DCUN1D1 in 37 %, FGFR1 in 19 % whereas none of the cases were amplified in ERBB2 gene. The majority of the samples exhibited amplification in at least one gene while half of them displayed concurrent amplification of two/three genes. Interestingly, 93 % of the FGFR1 amplified cases were also found co amplified with DCUN1D1 and/or BCL9 genes. Linear correlations were found between BCL9 and DCUN1D1 as well as BCL9 and FGFR1 gene amplification. BCL9 and DCUN1D1 genes' amplification was correlated with poorly differentiated tumors (p = 0.035 and p = 0.056 respectively), implying their possible role in tumor aggressiveness. This is the first study, to the best of our knowledge that examines the correlation of DCUN1D1 and BCL9 genes relative copy number amplification with molecular alterations and clinicopathologic characteristics of squamous cell lung cancer tissue samples. Our findings show concurrent amplification of genes in different chromosomes, with possible involvement in tumor aggressiveness. These results support the complexity of LUSC tumorigenesis and imply the necessity of multiple biomarkers / targets for a more effective therapeutic result in LUSC.

Estimating Copy-Number Proportions: The Comeback of Sanger Sequencing.

Determination of the relative copy numbers of mixed molecular species in nucleic acid samples is often the objective of biological experiments, including Single-Nucleotide Polymorphism (SNP), indel and gene copy-number characterization, and quantification of CRISPR-Cas9 base editing, cytosine methylation, and RNA editing. Standard dye-terminator chromatograms are a widely accessible, cost-effective information source from which copy-number proportions can be inferred. However, the rate of incorporation of dye terminators is dependent on the dye type, the adjacent sequence string, and the secondary structure of the sequenced strand. These variable rates complicate inferences and have driven scientists to resort to complex and costly quantification methods. Because these complex methods introduce their own biases, researchers are rethinking whether rectifying distortions in sequencing trace files and using direct sequencing for quantification will enable comparable accurate assessment. Indeed, recent developments in software tools (e.g., TIDE, ICE, EditR, BEEP and BEAT) indicate that quantification based on direct Sanger sequencing is gaining in scientific acceptance. This commentary reviews the common obstacles in quantification and the latest insights and developments relevant to estimating copy-number proportions based on direct Sanger sequencing, concluding that bidirectional sequencing and sophisticated base calling are the keys to identifying and avoiding sequence distortions.

Aberrant expression of thyroid transcription factor-1 in meningeal solitary fibrous tumor/hemangiopericytoma.

Meningeal solitary fibrous tumor (SFT) and hemangiopericytoma (HPC) were categorized as the same entity in the World Health Organization (WHO) 2016 classification of tumors of the central nervous system (CNS). Although NAB2-STAT6 fusion protein can be used to distinguish most of SFT/HPC from the other sarcomas, additional biomarkers were requested to separate meningeal SFT/HPC from meningioma and the molecular pathological difference between meningeal SFT/HPC and extra-CNS SFT/HPC still remains unclear. In this study, we evaluated the expression of TTF-1 in 67 meningeal SFT/HPC, 62 extra-CNS SFT/HPC and 201 meningiomas samples with immunohistochemistry (IHC) assays. The results showed that TTF-1 was detected in 23 of 67 (34.3%) meningeal SFT/HPC, 3 retroperitoneum SFT/HPC and none of meningiomas. Meanwhile, the copy number variation and mRNA expression of TTF-1 were measured by real-time quantitative PCR (qPCR) in meningeal SFT/HPC. These results demonstrated that TTF-1 protein expression level was significantly correlated with its transcription level, but independently related to the gene copy number variant. In conclusion, our study suggested that a large proportion of meningeal SFT/HPC was positive to TTF-1, while very few extra CNS SFT/HPC cases and no meningiomas were stained. So TTF-1 has value as an auxiliary diagnostic marker for meningeal SFT/HPC.

Consequences of aneuploidy in human fibroblasts with trisomy 21.

An extra copy of chromosome 21 causes Down syndrome, the most common genetic disease in humans. The mechanisms contributing to aneuploidy-related pathologies in this syndrome, independent of the identity of the triplicated genes, are not well defined. To characterize aneuploidy-driven phenotypes in trisomy 21 cells, we performed global transcriptome, proteome, and phenotypic analyses of primary human fibroblasts from individuals with Patau (trisomy 13), Edwards (trisomy 18), or Down syndromes. On average, mRNA and protein levels were increased by 1.5-fold in all trisomies, with a subset of proteins enriched for subunits of macromolecular complexes showing signs of posttranscriptional regulation. These results support the lack of evidence for widespread dosage compensation or dysregulation of chromosomal domains in human autosomes. Furthermore, we show that several aneuploidy-associated phenotypes are present in trisomy 21 cells, including lower viability and increased dependency on serine-driven lipid synthesis. Our studies establish a critical role of aneuploidy, independent of triplicated gene identity, in driving cellular defects associated with trisomy 21.

A statistical framework for determination of minimal plasmid copy number required for transgene expression in mammalian cells.

Plasmid DNA (pDNA) has been widely used for non-viral gene delivery. After pDNA molecules enter a mammalian cell, they may be trapped in subcellular structures or degraded by nucleases. Only a fraction of them can function as templates for transcription in the nucleus. Thus, an important question is, what is the minimal amount of pDNA molecules that need to be delivered into a cell for transgene expression? At present, it is technically a challenge to experimentally answer the question. To this end, we developed a statistical framework to establish the relationship between two experimentally quantifiable factors - average copy number of pDNA per cell among a group of cells after transfection and percent of the cells with transgene expression. The framework was applied to the analysis of electrotransfection under different experimental conditions in vitro. We experimentally varied the average copy number per cell and the electrotransfection efficiency through changes in extracellular pDNA dose, electric field strength, and pulse number. The experimental data could be explained or predicted quantitatively by the statistical framework. Based on the data and the framework, we could predict that the minimal number of pDNA molecules in the nucleus for transgene expression was on the order of 10. Although the prediction was dependent on the cell and experimental conditions used in the study, the framework may be generally applied to analysis of non-viral gene delivery.